Lya Escape from z ~ 0.03 Star-Forming Galaxies: the Dominant Role of 

Outflows 



Aida WofFord^, Claus Leitherer^ and John Salzer^ 
^ Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218; 

wof f ordSstsci . edu 

"^Astronomy Department, Indiana University, Swain West 408, 727 East Third Street, 

Bloomington, IN 47405 

ABSTRACT 

The usefulness of H i Lya photons for characterizing star formation in the distant 
universe is Hmited by our understanding of the astrophysical processes that regulate 
their escape from galaxies. These processes can only be observed in detail out to a 
few xlOOMpc. Past nearby {z < 0.3) spectroscopic studies are based on small samples 
and/or kinematically unresolved data. Taking advantage of the high sensitivity of HST^s 
COS, we observed the Lya lines of 20 Ha-selected galaxies located at {z) = 0.03. 
The galaxies cover a broad range of luminosity, oxygen abundance, and reddening. 
In this paper, we characterize the observed Lya lines and establish correlations with 
fundamental galaxy properties. We find seven emitters. These host young (< lOMyr) 
stellar populations, have rest-frame equivalent widths in the range 1 — 12 A, and have 
Lya escape fractions within the COS aperture in the range 1 — 12 %. One emitter has a 
double-peaked Lya with peaks 370 km s~^ apart and a stronger blue peak. Excluding 
this object, the emitters have Lya and Oi A1302 offsets from Ha in agreement with 
expanding shell models and LBG observations. The absorbers have offsets that are 
almost consistent with a static medium. We find no one-to-one correspondence between 
Lya emission and age, metallicity, or reddening. Thus, we confirm that Lya is enhanced 
by outflows and is regulated by the dust and Hi column density surrounding the hot 
stars. 

Subject headings: galaxies: evolution — galaxies: ISM — galaxies: starburst — galaxies: 
stellar content — ultraviolet: galaxies 



1. INTRODUCTION 



T he potential of Hi Lya to i dentify young galaxies at high redshift was noted over 40 years 
ago bv [Partridge &: Peebles! (|l967l ). Lya corresponds to the n = 1 — 2 transition and is in principle 
the strongest emission line originating in an Hii region. In the absence of dust, the Lya /Ha 
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line intensity ratio is expec ted to be in the range 7 — 12 for recombination theory cases B-A 
([Osterbrock &: Ferlandl l2006l ) . and the Lya equivalent width, EW(Lya), is expected to be in the 
range 100 — 300 A for metallicities between Z = 0.020 — 0.004, a nor mal initial mass functi on, a 
constant star formation rate of 1 M0/yr~^, and ages below 10 Myr ( Verhamme et al.ll2008l ). In 
the case of a single burst of star formation, the range of EW(Lya) is ~ — 300 under the same 
conditions. A few times 10^ — 10^ yr af t er a burst of star form ation, EW(Lya) drops to zero 



(jCharlot fc Falll Il993l : IValls-GabaudI Il993l : IVerhamme et al.l 120081 ) . The above predictions do not 
account for the underlying stellar contribution at Lya. This contribution is not a concern in galaxies 
dominated by populations of O and early B stars. For re f erenc e, the mean Lya equivalent width 
of the eight early and mid-0 supergiants in lBouret et al.l ( 20121 ) is ~1 A. On the other hand, the 
correction can be significant for stellar populations dominated by A and late B stars. In addition, 
the above predictions do not account for the effects of Lya radiation transfer in the medium 
surrounding the starburst. Due to its resonant nature, a Lya photon is absorbed and re-emitted 
multiple times in H i. Resonant trapping reduces the mean free path of the Lya photons increasing 
significantly their probability of being destroyed by dust, shifted in frequency, or transformed to 
two-photon emission. On the other hand, some effects are favorable to the escape of Lya photons 
from galaxies. Neutral gas outflows Doppler-shift the Lya photons from the optically t hick line core 



to the optically thin line wi ngs making it possible for some Lya photons to escape (iKunth et al 



19981 : 1 Verhamme et al.l l2006l). Lya photons may also escape thro ugh holes of sufficiently low dust 
and Hi column densities (jGiavalisco et al.lll996l : lAtek et al.ll2009l ). 



In spite of the various Lya attenuation effects, today, star-forming galaxy candidates have been 
detected out to redshift z ~ 10 using various techniques. Galaxies with redshifts of z ~ 2 — 5 are 
detected via the drop-out technique base d on the opacity at the Lyman- limit of the galaxy and the 
intervening interg alactic medium (IGM; iGiavaliscd l2002l : ISteidel et al.l 120031 : IShapley et al.l 120031 : 
Ouchi et al.ll2004l ). Gal axies with z > 6 are detected via the drop-out technique based on the opac- 
ity of the IGM at Lya (|Bouwens et al.1 boO?! : loesch et al.l[201ol : iBouwens et al.1 [2011al lbll . Galaxies 
with z ^2 — 7 are also detected using deep i iarrow band surveys that target the Lya enr i ssion- 



line (Cowie &: Hu 



1998 



Hu et al 



1998 



2004 



arrow bang surveys tnat target tne .Lya en u ssion- 
2010'; 'Rh oads et all bo od: 'Ou chi et ahl boo s'. '2008"; 



Taniguchi et al.lbo05l : IVenemans et al.lbo05l : lNilsson.2007i : INilsson et al.i,2009, : iGronwall et al..,2007l 



Finally, galaxies at the highest redshifts are detected using gravitationally lensin g clusters (lOta et al 



2OI2II. Hundreds of galaxy candidates have been spectroscopically confirmed (ISteide. 



Rhoads et al. 



2003 



Dawson et al.ll2007l : llye et al.ll2006l : iCowie et al.ll2010l : iBradac et al. 



1996 



few galaxies show impressive l y high rest-frame equiva lent widths of EW(Lya) >150 A ( Malhotra &: Rhoads 



2002 



Shimasaku et al 



2006: 



ence of an AGN dcharlot fc FalL,199.C 



et al ■ 

2012l'l. and a 



Finkelst ein et al.ll2008l ). which could alternatively be due to the pres- 



In principle, after removing the observational and astr qphysical biases, Lya can be used for var- 



ious applications. 1) To probe cosmic sta r formation rates (jKudritzki et al 



Gronwall et al 



2007 



200c 



Fujita et al. 



2003 



Zheng et al.ll2012l ). 2) To probe the ioniza.tion fraction of the intergalactic 



medium during the final stages of re-ionization ([Kashikawa et al 



2006 



Malhotra &: Rhoads 



2006 
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Dawson et a. 



Ouchi et al 



200^ 



200^ 



2001 



200 



2004 l2005l : 



2007 : IZheng et al.ll20ld: iBradac et al 



Venemans et al 



20121') . 3) To probe the large-scale structure 



20021 . l2004l : IWang et alJl2005l : IStiavelli et al 



Gawiser et al.ll2007l) . 4) To identify po tential hosts of population III star formation (jMalhotra Sz Rhoads 



redshift galaxie s 



Shimasaku et alJbood: ISch aerei 2007 ). 5) And to c onstrain the nature and e volution o f high 



Giavahsco ( 200^1 : iNilsson et all j^OoJ); IPinkelstein et alJ ^20od ): Ishaplevl (|201lh : 



Malhotra et alJ (l2012l ): iGonzalez et alJ (l2012l ). However, the use of Lya as a cosmological tool is 
limited by our understanding of the astrophysical processes that regulate the escape of Lya photons 
from galaxies and affect them on their way to telescopes. 

The Lya escape problem can be approached from two complementary angles, doing statistically 
significant studies of the global properties of the sources and their cosmological context, or studying 
individual sources in detail. The former is better done in the distant universe, whereas the latter 
can only be accomplished locally, where the effect of the intergalactic medium is a non-issue, the 
kinematical and spatial resolutions are higher, and ancillary data is abundant. The challenges of 
studying Lya locally are the necessity of a UV space-based observatory optimized for studying 
Lya; contamination with the geocoronal Lya line, which saturates the spectra at very low redshift; 
and contamination with the Milky Way Galaxy Lya-line trough. Indeed, the interstellar H I in our 
Galaxy produces a damped absorption along many sight-lines, which hides the potential Lya line 
from star- forming galaxies with redshifts below a few xlOO km s~^. 



There are several observational Lya studies of galaxies at low-redshift (z<0.3). iGiavalisco et al 



(jl996l ) studied low resolution spectra of 21 galaxies taken with the International Ultraviolet Ex- 
plorer (lUE) and studied the effects of reddening and metallicity on EW(Lya), and Lya/H/3. They 
found no correlation between these quantities and the UV extinction, the Balmer decrement, or 
the oxygen abundance. They concluded that rather than the amount of dust, the ISM geometry 
determines the Lya escape fraction. Unfortunately, t he spectral resoluti on of lUE is insufficient 
to study the effect of gas outflows on the Lya escape. iKunth et al.l (|l998l ) studied high-dispersion 
Hubble Space Telescope (HST) GHRS spectra of eight nearby Hll galaxies. They found four emit- 
ters, all of which have a velocity offset between the Hi and the Hil gas of up to 200 km s~^, while 
the non-emitters sho w no velocity offset. However, their study suffers from small number statistics. 
Scarlata et al.l (120091 ) analyzed the optical spectra of a sample of 31 {z) ~ 0.3 Lya emitters iden- 
tified by Deharveng et al. (2008) in low resolution (~ 8 A) Galaxy Evolution Explorer ( GALEX) 
data. They tried to reproduce the stellar-absorption corrected Lya/Ha and Ha/H/3 ratios using 
different dust geometries and assumptions for the Lya scattering in Hi. They found that the ratios 
are well reproduced by a clumpy dust distribution, while a uniform dust screen model results in 
Lya/Ha larger than the case B recombination theory value. Unfortunately, the spectral resolution 
of GALEX is insufficient for kinematical studies, and the latter a uthor s were unable to assess the 
relative importance of dust geometry and outfiows. lAtek et al.l (|2009l ) used a 3D Lya radiation 
transfer code to model HST observations of the very metal-poor dwarf galaxy I Zw 18, which shows 
Lya in absorption. Their analysis shows that it is possible to transform a strong Lya emission of 
into a damped absorption, even with low extinction, if the Hi column density 
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is large. Finally, lOstlin et alj (|2009l ) used HST to produce Lya, Ha, and UV continuum maps 
of six galaxies covering a wide range in luminosity and metallicity, including known Lya emitters 
and non-emitters. They found that the bulk of Lya emerges in a diffuse component resulting from 
scattering events. They also found the simultaneous presence of Lya in emission and absorption 
within spatial scales ranging from a few xlO pc to a few xlOO pc. 

Although these studies have addressed main factors affecting the Lya escape, we still do not 
understand the relative importance of starburst phase, gas kinematics, gas geometry, and dust 
extinction, in determining the profile and strength of the Lya line from star-forming galaxies. We 
took advantage of the high sensitivity and medium resolution of the Cosmic Origins Spectrograph 
(COS) aboard HST in order to obtain Lya spectroscopy of a sample of 20 galaxies located a mean 
redshift of 0.03 and covering a broad range of luminosity, oxygen abundance, and reddening. We 
used these and ancillary data for establishing correlations between Lya and fundamental galaxy 
properties. In section 2 we describe the sample, the observations, and the ancillary data; in section 
3 we present our analysis; in section 4 we discuss our results; and in section 5 we summarize and 
conclude. 



2. SAMPLE, OBSERVATIONS, AND ANCILLARY DATA 

2.1. Sample 

The galaxies in our sanipl e were Ha-selected from the Kitt Peak International Spectroscopic 
Survey data release (KISSbFI. Isalzer et al.l [2OO1I '). The sample is composed of 12 irregular and 



eight spiral galaxies distributed in the redshift range z = 0.02 — 0.06. At these redshifts, we avoid 
contamination with the geocoronal and Milky Way (MW) Lya lines. We cover a factor of three 
in redshift, which makes it possible to study the effect of distance on the observed properties of 
Lya. In addition, we cover a broad range in oxygen abundance, 12 -|- log (O/H) = 7.9 — 9.1, which 
makes it possible to study the effect of metallid;|l on the Lya escape. Metallicity is expected to 
correlate with the dust content, and dust destroys Lya photons. We also cover a wide range in 
Ha equivalent width, EW{Ka) = 27 — 578 A. The value of EW{ Y{q) is sensitive t o the starburst 



evolutionary phase and to the initial mass functior i (IMF), e.g., seelLeitherer et al.l (|l999l ). Finally, 



we cover over an order of magnitude in GALEX ([Martin et al.ll2005l ) FUV continuum luminosity. 



log L(1500) = 39.0—40.2 erg s~^ (uncorrected for reddening). For a given star formation history 
(SFH), age, and IMF, L(1500) is proportional to the stellar mass in the population dominating 
the FUV emission of the galaxy. Therefore, -L(1500) can be used to determine the star formation 
rate of this population (SFR = stellar mass / time since beginning of star formation). Table [T] 
summarizes the properties of the sample. The morphological types are based on SDSS images. The 



^The R in "KISSR" stands for data with red spectra, i.e., from 6400 - 7200 A 

^In this paper, metaUicity is the mass ratio of aU elements excluding H and He to all elements. 
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Ha equivalent widths are from the SDSS spectra when available or from the KISS data otherwise, 
as indicated in column (10) of Table [TJ 



2.2. Spectral Class of the Galaxies 



We are interested in galaxies where the principal photoionization source is the strong UV 
radiation from OB stars. This helps in the interpretation of the observed Lya properties because 
we do not need to worry about non-stellar ionizing sourc es such as active galactic nuclei (AGNs ), 
which account for some of the high Lya EWs at z > 2 (ICharlot &: Fall! Il993l : IZheng et al.l l20ld ) , 
and low-ioniza tion nuclear emission-line regions (LINERs), whose nature is not well understood 
(jTanakall2012l ). We verified the spectral class of our galaxies using the so-called BPT diagnostic 
diagram (Baldwin, Phillips & Terlevich 1981). This was not possible for KISSR 271, for which 
an optical spectrum was unavailable. As Figure [T] shows, most galaxies are dominated by star 
formation, but KISSR 1084 may be of composite type. 



2.3. Observations 



Our observations are part of COS-GTO programs 11 522 and 12027 (P . I. J. G reen). Descriptions 
of the COS and its on-orbit performance can be found in lOsterman et al.l (j201ll ). We used one orbit 
with HST to observe the brightest FUV knot detected by GALEX in each galaxy. We acquired 
the targets using the ACQ/IMAGE mode and mirror A. Therefore, we obtained two images of 
each galaxy in the observed-frame wavelength range 1700 — 3200 A. One after the initial telescope 
pointing, and the second after the final telescope pointing. We then used the G130M grating to 
obtain a spectrum of each target in the observed-frame wavelength range 1138-1457 A. We achieved 
continuous wavelength coverage by using central wavelengths, A1291 and A1318 at the default focal 
plane offset position (FP-P0S=3). Table [2] lists the epochs and exposure times corresponding to 
these central wavelengths. The spectra were taken through the circular Primary Science Aperture 
(PSA), which is 2.5" in diameter. They contain the aberrated light from objects located up to 2" 
or up to 0.9 — 2.4 kpc from the aperture's center, depending on the distance to the galaxy. The 
data were pr ocessed with CALCOS v2.13.6 and combined with the custom IDL co-addition routine 
described in iDanforth et al.l (|2010l ) . 



2.4. COS Images 

Figure [2] shows the NUV target acquisition (TA) confirmation images of the galaxies with the 
COS aperture overlaid. The center of the aperture is indicated by the small black cross. In 12 
cases, a clear source centroid associated with young star clusters is present. However, in eight 
cases, there is no clear centroid of the emission. This is due to the short exposure times of the TA 
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Table 1. Sample Properties 



Galaxy 


Morph. 


z 


d 


Mb 


E(B-V)g 


12-t-log(0/H) 


EW(Ha) 


logLisoo 


Note 








Mpc 


mag 


mag 


dcx 


A 


erg/s/A 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


40 


Irr 


0.026 


112 


-18.0 


0.03 


8.18 


58 


39.29 


1 


J.UO 


Irr 


0.024 


101 


-17.6 


0.01 


8.17 


108 


39.22 


1 
± 


178+ 


Irr 


0.057 


246 


-20.1 


0.01 


8.52 


100 


40.13 


1 


182 


Irr 


0.022 


96 


-17.8 


0.01 


8.34 


114 


39.33 


1 


218 


S 


0.021 


90 


-19.4 


0.01 


8.87 


66 


39.65 


1 


242* 


Irr 


0.038 


162 


-19.4 


0.01 


8.38 


435 


40.25 


1 


271* 


Irr 


0.023 


97 


-17.6 


0.02 


8.39 


29 


39.13 


2 


298* 


S 


0.049 


210 


-20.0 


0.02 


9.07 


51 


40.06 


1 


326* 


s 


0.028 


120 


-20.0 


0.02 


8.74 


78 


39.65 


2 


1084* 


s 


0.032 


137 


-20.5 


0.05 


8.79 


87 


39.7 


1 


1567* 


s 


0.045 


194 


-20.0 


0.02 


8.75 


29 


40.18 


2 


1578* 


Irr 


0.028 


120 


-19.6 


0.02 


8.14 


576 


40.29 


1 


1637 


Irr 


0.035 


149 


-19.6 


0.01 


8.55 


34 


39.92 


1 


1785- 


Irr 


0.021 


90 


-17.6 


0.01 


7.99 


199 


39.07 


1 


1942 


Irr 


0.039 


169 


-19.3 


0.03 


8.59 


25 


39.83 


1 


2019 


Irr 


0.034 


146 


-17.9 


0.02 


7.96 


465 


39.48 


1 


2021 


S 


0.04 


171 


-19.6 


0.02 


8.79 


79 


39.93 


1 


2023 


s 


0.036 


155 


-19.4 


0.02 


9.00 


58 


39.77 


1 


2110 


Irr 


0.025 


108 


-17.3 


0.01 


7.86 


511 


39.36 


2 


2125 


S 


0.025 


108 


-18.7 


0.01 


8.76 


85 


39.31 


1 


Min 




0.021 


90 


-20.5 


0.01 


7.86 


25 


39.07 




Max 




0.057 


246 


-17.3 


0.05 


9.07 


576 


40.29 




Note. 


— (1) KISSR ID. We mark the Ly, 


3: emitters with asterisks and the furthest/nearest 


galaxies 



with plus/minus signs. The last two rows give the minimum and maximum values in each column. (2) 
Morphological type. S=spiral. Irr=irregular. (3) Redshift of Ha line. (4) Distance computed from the 
redshift using d = cz/Hq, where c is the speed of light, z is the redshift, and we adopt [^;^, Qmi ^^o] = [0.7, 
0.3, 70 k m s"-*^ Mpc~^1 . (5) A bsolute blue magnitude. (6) Color excess b ased on the Milky W ay extinction 
maps of ISchlegel et aU lll998l) and th e redd ening law with Rv = 3.1 of iFitzpatricld l ll999l ). (7) Oxygen 
abundance derived as in ISalzer et al. I 1I2OO5I ). (8) Ha equivalent width. (9) Luminosity at 1500 A from 
GALEX, uncorrected for reddening. (10) Source of the ancillary optical spectrum. 1=SDSS. 2=Same as 
in Salzer et al. 2005. 
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Fig. 1. — Baldwin, Phillips & Terlevich diagnostic diagram of the spiral (filled symbo l s) an d 
irregular (open symbols) galaxies in our sample. We adopt the scheme of iKewley et alJ (|2006l ). 
where the principal excitation mechanism in the galaxy is star formation if the galaxy lies below 
the dotted curve, an AGN if the galaxy lies above the solid curve, and a mixture of star formation 
and Seyfert nucleus, or star formation and LINER, if th e galaxy lies betw e en th e dotted and solid 
curves. The dotted curve cor responds to e q uation 1 in iKauffmann et al.l (j2003l ). The solid curve 
corresponds to equation 5 in iKewlev et al.l (j200ll ) . The symbols of the Lya emitters are labeled 
with the KISSR ID of the galaxy. 
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images and in two cases, i.e., for KISSR galaxies 40 and 1567, to TA failures. In spite of the TA 
failures, spectra were taken. The galaxies that turned out to have net Lya emission are marked 
with asterisks. Hereafter, we will call the latter galaxies the emitters. Note that emitters KISSR 
271, 326, and 1567, do not show source centroids. Also note that in spit e the TA failure, K ISSR 
1567 turned out to be an emitter. This is in agreement with the finding of lOstlin et al.l (j2009l ) that 
the Lya emission is not necessarily coincident with regions of high FUV surface brightness. 



2.5. Ancillary Images 

We obtained information about the morphology and orientation of the galaxies from archival 
images. At the time of our analysis, there were no HST images available for our targets. The 
galaxies were imaged by GALEX in the FUV at low spatial resolution (~ 0.5"). Unfortunately, in 
these li nages, the galaxie s look like large blobs. However, we found useful Sloan Digital Sky Survey 



(SDSS, lYork et al.l l200d ) images of all the targets in our sample. These are shown in Figure [3l 
where we have overlaid the COS footprint at the approximate location of the telescope pointing 
and marked the emitters with asterisks. 

Figure [3] shows that most spiral galaxies in our sample have low inclinations except KISSR 
2125. While some low-inclination spirals have Lya in emission, others don't. KISSR 1084, the 
spiral of composite spectral class is an emitter, and as expected, the highly inclined spiral is a 
non-emitter. Our two most metal-poor irregulars, KISSR 2019 and 2110, are Lya absorbers, but 
irregulars KISSR 242 and 1578 are Lya emitters. In summary, there are emitters and absorbers 
among both morphological types. Finally, note that the COS aperture only covers an area of 
0.9 — 2.4 kpc in radius, depending on the distance to the galaxy. This needs to be considered when 
comparing with Lya observations at higher redshifts , since at higher redshift, more of the galaxy is 
enclosed within the aperture and Ostlin et al. (j2009l ) found that the bulk of the Lya photons from 



nearby galaxies are in the diffuse halo. 



2.6. Ancillary Spectra 



We used ancillary optical spectroscopy of the galaxies for determining their redshift, spectral 
class, metallicity, reddening, SFR, a nd Wolf-Rayet (WR) s tar content. Sixteen targets have spectra 
from the SDSS seventh data release ([Abazajian et al.ll2009l ). The SDSS spectra are of higher spectral 
resolution than the KISSR spectra, and they were taken through a circular fiber of 3" in diameter 
that closely matches the size of the COS aperture. The difference between the SDSS and the COS 
pointings is relevant to the interpretation of the Lya/Ha ratio, and to the comparison of the strength 
of Lya and the value of EW{Ha). Indeed, the COS aperture encompasses only a portion of each 
galaxy. Table [3] lists the KISSR identifier (Column 1, the emitters are marked with asterisks), the 
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40 UK 17S 3S 318 




Fig. 2.— Near-ultraviolet (1700 - 3200 A) COS TA confirmation images. North is up and east is 
to the left. The numbers are the galaxy identifiers. The circles represent the 2.5" in diameter COS 
PSA aperture. The emitters are marked with asterisks. In the rainbow scale, red represents the 
region of maximum counts, while pink represents regions of low counts. 
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K1SSH_1C« IKISSR i7B. lK|3aH_1K. K]Sai_21$. 



<iSf5HJ4? Ug3sn_Z7i jqssH^aH. 

ji^i^AJH*ms^7 [jmM&re^»aisi.i J4s»4«i7»»»«4u 



o 



« 






6 




WS3Ft_t7S5 Kissa_ws. 



J134a44.4f*4f1MB.S 



® 





Fig. 3. — Sloan Digital Sky Survey 25"x25" u-, g-, r-, i-, and z-band composite images of the 
galaxies. The KISSR and SDSS galaxy identifiers are provided. North is up and east is to the 
left. The Lya emitters are marked with asterisks. The yellow circles show the COS footprint and 
pointing. 
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right ascension (RA, J2000) and declination (DEC, J2000) of the first HST pointing (Columns 2 
and 3), the RA and DEC of the SDSS spectroscopic pointing (Columns 4 and 5), and the difference 
between the COS and the SDSS pointings (Columns 6-8). As Table [3] shows, on average, the SDSS 
pointings are within 0.3" of the COS initial pointings, with a standard deviation of 0.4". Since the 
diameter of the COS aperture is 2.5", we do not expect the difference in pointings to constitute a 
major issue. For KISSR galaxies 326, 1567, and 2110, which do not h ave SDSS spectra, we were 
able to recover the optical 2" -wide slit spectra of I Jan gr en et al.l ()2005l ). Unfortunately, we do not 
have an optical spectrum for KISSR galaxy 271. 



3. ANALYSIS 

3.1. Wavelength Zero Point and EfFective Spectral Resolution 

The COS spectra include the intrinsic lines Lya and a number of low-ionization (LIS) and 
higher-ionization interstellar lines. We are interested in comparing the offsets of these lines with 
respect to Ha. This is for two reasons. First, as mentioned in the introduction, velocity offsets 
between the Hll gas traced by Ha and the Hi gas traced by the LIS lines have been showed to 
favor the escape of Lya from galaxies. Second, in expanding-shell radiation tr ansfer models the 



Lyg emission line is redshifted to a few times the shell expansion velocity (e.g., IVerhamme et al 



20061 ). A similar offset is observed at high redshift when comparing the Lya redshift to the blueshift 



of the LIS lines, which give the expansion velocity (e.g. IShapley et al.ll2003l ). Thus, offset measure- 
ments are useful for establishing how the Lya strength relates to the interstellar gas kinematics 
and for constraining radiation transfer models. The accuracy of the velocity offset measurements 
depends on the accuracy of the wavelength zero point and on the effective spectral resolution of 
the data. 

In ultraviolet spectroscopic studies of low-redshift galaxies, it is common to use geocoronal 
emission lines to check if the wavelength zero point is shifted due to instrumental effects, and Milky- 
Way metal absorption lines to check if the wavelength ze ro point is shifted due to the asymmetry 



of the target in the aperture (e.g., see iKriss et al.ll201ll ). Figures SHS] show the COS spectra of 



the galaxies. They are contaminated with some or all of the geocoronal emission lines H i A1216, 
N1AI2OO, and/or Oi A1302, 1305; and with the Milky Way (MW) absorption lines. Sin A1190, 
1193, 1260, 1304, Sim A1206, Oi 1302, 1305, and/or C 11 A1334, 1336. We measured the centroids 
of the contaminating lines using a custom routine developed by COS science team member K. 
France. The routine takes into consideration the COS line spread function. Table |4] lists the mean 
wavelength offsets of the geocoronal and MW lines. We find that the mean velocity offsets of the 
geocoronal and MW lines are in the ranges —33 to 18 km s~^ and —1 to 110 km s~^, respectively. 



^The position of the final HST pointing is not recorded in the standard output from CALCOS v2.13.6. It is 
expected to be within 0.3" of the initial pointing based on results of blind pointing exposures. 
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Table 2. Epochs and Exposure Times of G130M Observations 



Galaxy Start Time Start Time Exp. Time Exp. Time 





1291 A 
(2) 






1318 A 
(3) 








\^s) loio A 


40 


2011/06/11 


02; 


;07; 


;49 


2011/06/11 


03; 


:26; 


:08 


1105 


1105 


108 


2011/07/31 


04; 


;47; 


;42 


2011/07/31 


05; 


:09; 


:24 


1105 


1105 


178 


2011/06/17 


20; 


;55; 


;13 


2011/06/17 


21 


;16; 


;55 


1105 


1105 


182 


2011/07/26 


08; 


;13; 


;51 


2011/07/26 


08 


;35; 


;33 


1105 


1105 


218 


2010/02/23 


07; 


;33; 


;10 


2010/02/23 


09; 


;09; 


;06 


1111 


1111 


242 


2009/12/26 


23; 


;56; 


;45 


2009/12/27 


00; 


;18; 


;33 


1111 


1111 


271 


2011/05/07 


17; 


;41; 


;45 


2011/05/07 


18; 


:03; 


:27 


1105 


1105 


298 


2010/01/01 


12; 


;41; 


;22 


2010/01/01 


13: 


;03; 


;04 


1105 


1105 


326 


2011/04/22 


13; 


;36; 


;31 


2011/04/22 


13; 


;58; 


;13 


1105 


1105 


1084 


2010/02/05 


05; 


;13; 


;19 


2010/02/05 


05; 


;35; 


;01 


1105 


1105 


1567 


2010/08/13 


05; 


;09; 


;19 


2010/08/13 


05; 


;31; 


;46 


1150 


1150 


1578 


2010/02/18 


22; 


;07; 


;27 


2010/02/18 


22; 


;29; 


;57 


1153 


1153 


1637 


2009/10/24 


09; 


;51; 


;17 


2009/10/24 


10; 


:13; 


:44 


1150 


1150 


1785 


2011/08/09 


10; 


;38; 


;26 


2011/08/09 


11; 


:00; 


:55 


1142 


1142 


1942 


2011/06/13 


22; 


;46; 


;13 


2011/06/13 


23; 


;08; 


:28 


1128 


1128 


2019 


2011/06/19 


15; 


;57; 


;20 


2011/06/19 


16; 


;19; 


;49 


1142 


1142 


2021 


2009/11/14 


00; 


;02; 


;53 


2009/11/14 


00; 


;25; 


;17 


1147 


1147 


2023 


2011/06/16 


14; 


;29; 


;53 


2011/06/16 


14; 


;52; 


;08 


1128 


1128 


2110 


2011/06/16 


16; 


;16; 


;21 


2011/06/16 


16; 


;38; 


:50 


1142 


1142 


2125 


2011/06/16 


17; 


;52; 


;32 


2011/06/16 


18; 


:15; 


:01 


1142 


1142 



Table 3. Comparison of COS and SDSS Pointings 



Galaxy 


COS RA 


COS DEC 


SDSS RA 


SDSS DEC 


ARA 


AD EC 


Atot 




deg 


deg 


deg 


deg 


deg 


deg 


arcsec 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


40 


185.59833 


29.44383 


185.59825 


29.44394 


8E-05 


-lE-04 


0.49 


108 


190.98042 


29.36964 


190.98044 


29.36961 


-2E-05 


3E-05 


0.13 


178 


195.42306 


29.38118 


195.42304 


29.38135 


2E-05 


-2E-04 


0.62 


182 


195.60688 


28.85808 


195.60688 


28.85807 


-5E-06 


lE-05 


0.05 


218 


197.31725 


29.36739 


197.31703 


29.36766 


2E-04 


-3E-04 


1.24 


242* 


199.01625 


29.38161 


199.01632 


29.38169 


-7E-05 


-8E-05 


0.38 


271* 


200.42000 


28.88308 












298* 


202.45746 


29.57972 


202.45748 


29.57972 


-2E-05 


2E-06 


0.08 


326* 


204.81708 


28.87361 












1084* 


252.27192 


29.75878 


252.27195 


29.75878 


-3E-05 


-2E-06 


0.12 


1567* 


201.22750 


43.79889 












1578* 


202.18354 


43.93069 


202.18357 


43.93070 


-3E-05 


-6E-06 


0.10 


1637 


205.50002 


42.76004 


205.50038 


42.76000 


-4E-04 


4E-05 


1.30 


1785 


216.62029 


43.47083 


216.62029 


43.47082 


2E-06 


lE-05 


0.05 


1942 


230.70200 


43.77347 


230.70197 


43.77347 


3E-05 


2E-06 


0.11 


2019 


236.43554 


44.26439 


236.43552 


44.26439 


2E-05 


-lE-06 


0.08 


2021 


236.46992 


44.26322 


236.46991 


44.26322 


7E-06 


2E-06 


0.03 


2023 


236.67746 


43.24367 


236.6774 


43.24359 


6E-05 


8E-05 


0.35 


2110 


242.04625 


43.63175 












2125 


242.58508 


43.00975 


242.58507 


43.00975 


lE-05 


OE+00 


0.05 
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The MW lines arise from absorption in halo clouds that may have non-negligible heliocentric 
velocities. For this reason, we analyzed the Hi 21 cm emission from intermediate velocity clouds 



(IVC s) along the lines of sight to our targets. The data are from the LAB H i survey (jKalberla et al 



20051 ). In Figure [9] we compare the centroids of the Hi IVC emission lines and the MW Sin A1990 
absorptions. The Siii absorptions are corrected for the wavelength offset given by the geocoronal 
lines. The Hi 21 cm emission is the mean within the radio beam at each velocity. We chose 
Sin A1990 because it is the only MW line that is not blended with other lines for any of the 
targets. The Hi component near km s~^ originates in the MW disk. We find that the maximum 
IVC cloud velocity is ~ —50 km s~^ in the local standard of rest frame (e.g., KISSR 40, 1567, 
and 1578). Although the comparison between the 21 cm emission and the COS data is difficult 
in cases where the signal to noise of the COS data is low, in general, the 21 cm emission and the 
Sin absorption overlap. However, in some of the highest signal-to-noise cases (e.g., KISSR 242, 
1578, 2019, and 2110), the minimum of the absorption is blueward of the bluest 21 cm emission 
peak. Note that the radio beam is large in comparison with the COS aperture, which makes this 
wavelength calibration method imperfect. Since large offsets between the 21 cm and the COS data 
are unlikely to be due to shifts in the wavelength zero point, we decided to only correct the COS 
data for the small shift in the geocoronal lines. 

Regarding the spectral resolution, the more extended the source, the worse the effective spectral 
resolution, and the worse the accuracy in the determination of line centroids. Figure [2] shows that 
some targets are more extended than others. We used the MW lines Sin A1190, 1193, 1260, Sim 
A1206, and/or Cii A1334 for determining the effective spectral resolution of the COS spectra. The 
dispersion of the G130M grating is 9.97 mA/pix and the resolution element is 6 pixels. Therefore, 
the nominal resolution is 0.06 A or ~16 km s~^. The effective spectral resolution of each spectrum 
is given by the FWHM of the MW absorption lines. For each galaxy. Column 4 of Table H] lists 
the mean effective spectral resolution. The mean FWHM of the MW lines is in the range 74 — 305 
km s~^. We can derive the centroids of the LIS lines with an accuracy of a tenth of the effective 
spectral resolution, as long as the signal to noise ratio is sufficient. 



3.2. Lya and ISM Kinematics 

Here we present the velocity offsets of the intrinsic lines mentioned in the previous section. 
Figure fTOl shows the Lya profiles of the galaxies. Seven objects have net Lya emission (emitters), 
eleven objects have damped Lya absorptions (absorbers), and two objects, KISSR 2021 and 2023 
have noisy Lya profiles blended with MW lines and that appear to be in absorption with perhaps 
some Lya emission. We measured the Lya centroids as follows. For the emitters, the centroid is 
the point of maximum Lya flux. For the absorbers, we used the center of the best-fit Voigt profile. 
The fits to the Lya absorptions were done by eye using a custom routine written by J. Tumlinson. 
We overlay the Voigt profiles in Figure [TOl 

Excluding KISSR 298, the emitters have P-Cygni like profiles with the emission component 
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Fig. 4.— COS spectra of galaxies KISSR 40 , 108, 178, and 182 from top to bottom (KISSR ID 
provided on the upper left of each panel). The spectra are corrected for the redshift and uncorrected 
for the reddening. They were first binned to 16 km s~^ (nominal resolution of the COS spectra) 
and then smoothed for clarity. We identify foreground (Foregr.) and intrinsic (Intr.) lines at the 
bottom of each panel. Geocoronal lines are marked with an encircled plus. The intrinsic lines are 
marked with vertical gray lines. We did not clip the geocoronal lines in order to show the extent 
to which they affect the profiles of nearby lines. 
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Fig. 5.— Similar to Figured but for galaxies KISSR 218, 242, 271, and 298 from top to bottom. 
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Fig. 6. — Similar to Figure|4]but for galaxies KISSR 326, 1084, 1567, and 1578 from top to bottom. 
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Fig. 7.— Similar to Figure H but for galaxies KISSR 1637, 1785, 1942, and 2019 from top to 
bottom. 
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Fig. 8.^ Similar to Figure H but for galaxies KISSR 2021, 2023, 2110, and 2125 from top to 
bottom. 
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Table 4. Wavelength Zero Point and Effective Spectral Resolution 



Galaxy 


/ \ 


/ \ 

\Vh,MW/ 


(i^ WHMmw/ 


[t WHMi) 




km s""*^ 


km 


km 


km 


(1) 


(2) 


(3) 


/'A\ 


(5) 


40 


-32 ± 21 


-78 ± 27 


179 ± 91 


287 


108 


-21 ± 5 


-1 ± 5 


109 ± 37 


351 


178 


-27 ± 5 


-62 ± 29 


197 ± 24 


332 


182 


-21 ± 13 


-30 ± 12 


132 ± 16 


210 


218 


10 


-26 ± 51 


173 ± 45 


391 


242* 


18 ± 5 


-34 ± 8 


110 ± 27 


227 


271* 


-17 ± 5 


-10 ± 43 


74 ± 16 


321 


298* 


15 


-49 ± 17 


75 ± 52 


337 


326* 


-33 


-90 ± 60 


107 


274 


1084* 


10 ± 2 


-55 ± 9 


185 ± 38 


446 


1567* 


-17 ± 2 








1578* 


-11 


-70 ± 4 


125 ± 18 


308 


1637 


2 


-69 ± 71 


305 ± 131 


359 


1785 


-14 ± 2 


-110 ± 4 


142 ± 10 


246 


1942 


-7 ± 21 


-95 ± 61 


179 ± 39 


345 


2019 


-16 ± 6 


-101 ± 13 


188 ± 30 


238 


2021 


-2 ± 4 


-56 ± 25 


174 ± 81 


293 


2023 


-15 ± 9 


-38 ± 72 


198 ± 94 


304 


2110 


-8 ± 15 


-77 ± 9 


138 ± 21 


186 


2125 


-7 ± 14 


-78 ± 14 


274 ± 74 


349 


Min 


-33 


-110 ± 4 


74 ± 16 


186 


Max 


18 ± 5 


-1 ± 5 


305 ± 131 


446 



Note. — (1) KISSR ID of the galaxy. We mark the emitters 
with an asterisk. The last two rows give the minimum and maxi- 
mum values in each column. In cases where more than one line was 
measured, we give the standard deviation of the various measure- 
ments, as the error. (2) & (3) Average heliocentric velocities, t'^, 
of geocoronal and MW lines, respectively. (4) &; (5) Mean FWHM 
of MW and intrinsic ISM absorption lines. Values of < 110 km 
s~^ in Column 5 are lower limits. 
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-200-100 100 200 -200-100 100 200 -200-100 100 200 -200-100 100 200 -200-100 100 200 

LSR Velocity (km s"^) 

Fig. 9. — IVC Hi 21 cm emission along the lines of sight to our targets (thick gray curves) and 
MW Sill A1190 absorptions (thin black curves). The spectral binnings are the nominal for the 
LAB and COS data, i.e., 1 km for the former and 16 km for the latter. We give the KISSR 
ID on the upper left of each panel. The seven galaxies with net Lja in emission are marked with 
asterisks. The x-axis gives the velocity in the local standard of rest. The y-axis gives the continuum 
subtracted flux divided by a constant. The two vertical dotted lines correspond to vi,sr = and 
50 km s~^. 
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redshifted to a mean velocity of {v{Lya)) = 172±49 km s'^. KISSR 298 is a face-on spiral with two 
Lya peaks separated by 370 km s~^, such that one peak is redshifted and the other is blueshifted and 
twice as strong. The absorbers have Lya troughs centered at a mean velocity of {v(Lya)) = 0ib38 
km s~^. The Lya mean velocities exclude KISSR 2021 and 2023 due to their noisy Lya profiles. 

The COS spectra include the Hi gas tracers Oi A1302, Sin A1190, and Cii A1334. They also 
include Si iii A1206, which traces gas in a higher ionization stage. Figure [TT] shows the profiles of 
the previous metal lines in cases where the lines are not severely blended with contaminating lines. 
We omit emitter KISSR 1567, because for this object the target acquisition failed and the only 
observable line is Lya. The figure shows that C ii A1334 and C ii A1336 are resolved in the highest 
signal to noise spectra, e.g., for KISSR 242 and 1578. 

For the emitters, the mean blueshifts of Oi A1302 and Cii A1334 are {v{Ol)) = -117 ± 47 
km s~^ and (f (CII)) = —94 ± 20 km s~^, respectively. For the absorbers, the corresponding mean 
velocity offsets are (v(OI)) = —23 it 28 km and (v(CII)) = —21 it 31, which is almost consistent 
with a static medium. Note that for spiral KISSR 2023, whose Lya profile is very noisy, O i A1302 
has a large blueshift (—153 km s~^), such that Lya emission would be expected. As previously 
mentioned, there appears to be some Lya emission in this galaxy which has poor signal to noise. 
Interestingly, the SDSS image of KISSR 2023 is very similar to that of double-peaked emitter 
KISSR 298. 

Our kinematical results are shown in Table \5\ which lists the KISSR ID (Column 1); the 
Lya profile shape (Column 2); the centroid of the Lya line (Column 3); the velocity offsets between 
Ha and intrinsic interstellar absorption lines Sill A1190, Sin A1193, Si iii A1206, Ol A1302, Sill 
A1304, and Cii A1334 (Columns 4-9); the mean velocity offsets and standard deviations excluding 
Si in A1206 (Columns 10 & 11); the difference between the velocity offset of Si III A1206 and the 
value in Column 10 (Column 12); and for the Lya absorbers, the Hi column density of the Voigt 
profile (Column 13). 

In summary, excluding the double emitter KISSR 298, we find that: i) Lya is in emission and 
redshifted when high velocity outflows corresponding to expansion velocities of ~ 100 km s~^ are 
present, ii) Lya is in absorption and at a velocity close to that of Ha when null or low velocity gas 
flows are present, and iii) the Lya redshift is higher than the blueshift of the LIS lines by a factor 
of 1.5 — 2. These results are further discussed in section 4. In addition, in cases were the centroid 
of Si III A1206 was measured, we flnd that this line tends to have a blueshift larger than that of 
the LIS lines by about 30 — 50 km s~^. 

3.3. EW{Lya), Lya escape fraction, and -L1500 

In order to characterize the strength of Lya, we computed the Lya luminosity L(Lya), equiv- 
alent width, and escape fraction within the COS aperture. Although the relevant quantity for 
cosmological studies is the global integrated escape fraction of the galaxy, our measurements are 



- 22 - 




Fig. 10. — Intrinsic Lya profiles of tlic galaxies corrected for redshift. We give the KISSR ID on the 
upper left of each panel. The seven galaxies with net Lya in emission are marked with asterisks. The 
spectra are binned to 16 km s"^. In addition, the spectra of the absorbers are smoothed for clarity. 
We overlay the Voigt profiles used for determining the Lya absorption centroids (gray curves). 
The vertical dotted lines mark the rest-frame positions of Lya. Note the double Lya emission and 
stronger blue peak in KISSR 298. Also note that the Lya profiles of KISSR 2021 and KISSR 2023 
are rather noisy. 
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Fig. 11. — Velocity offsets between metal absorption lines intrinsic to the galaxies and Ha. We 
include Oi A1302 (red), Sin A1190 (blue), Cii A1334 (green), and/or Si iii A1206 (gray), except 
when they are blended with contaminating lines. We give the KISSR ID on the upper left of each 
panel. The emitters are marked with asterisks. The spectra are binned to 16 km s~^ and smoothed 
for clarity. The vertical dotted lines mark a velocity offset of zero. For KISSR 298, which is a 
double-peaked Lya emitter, we mark the position of the point equidistant from the two Lya peaks 
with a vertical dashed line. 



Table 5. ISM Kinematics 



Galaxy 


Lya 


Lya 1216 


Si II 1190 


Si II 1193 


Si III 1206 


O 1 1302 


Si II 1304 


C II 1334 


(Av) 


(t{Av) 


1206- (Av) 


N(Hi) 




shape 


km 


km 


km 


km s^^ 


km 


km 


km 


km 


km 


km 


cm^'^ 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


40 


Ab 


-22 • • • 




-29 








-29 


-29 







5.0E-t-20 


108 


Ab 


• • • 




-57 


-96 


-48 




-50 


-52 


5 


44 


7.5E-I-20 


178 


Ab 


2 • • ■ 








-35 


-48 


-17 


-33 


16 




2.3E-I-21 


182 


Ab 


-37 ■ • • 




-46 


-47 


-34 


-33 


-47 


-40 


8 


7 


2.0E-I-21 


218 


Ab 


10 ■ ■ • 








-32 


-11 


23 


-7 


28 




9.0EH-20 


242 


P-Cygni 


170 ± 46 


-68 


-73 


-130 


-53 


-104 


-79 


-76 


19 


55 




271 


P-Cygni 


152 ± 20 








-177 




-113 


-145 


45 






298 


2 Em 


-118 ± 11, 
249 ± 17 






-185 


-162 


-192 


-100 


-151 


47 


33 




326 


P-Cygni 


251 ± 24 


-162 


-140 




-83 


-175 


-88 


-130 


42 






1084 


P-Cygni 


198 ± 20 


-169 


-86 


-146 


-109 


-125 


-78 


-113 


36 


33 




1567 


P-Cygni 


105 ± 41 






















1578 


P-Cygni 


158 ± 2 


-163 


-147 


-166 


-116 


-149 


-127 


-140 


19 


26 




1637 


Ab 


-50 • • • 


-68 


-65 


-60 


-54 


-69 


-9 


-53 


25 


6 


7.0E+20 


1785 


Ab 


50 ••• 








-24 


-61 


-42 


-42 


19 




l.OE-l-21 


1942 


Ab 


-50 ■ • • 


-7 


-2 




6 


-63 


-24 


-18 


27 




9.5E-I-20 


2019 


Ab 


50 ■ ■ • 


-54 


-64 


-37 


-25 




-59 


-50 


18 


-14 


3.5E-I-20 


2021 


Ab 


-50 


-38 


-28 


-19 


-8 


-18 


-14 


-21 


12 


-2 


2.5E-I-20 


2023 


Ab 


-50 








-153 


-62 


-171 


-128 


58 




l.OE-l-20 


2110 


Ab 


50 ••• 


-13 


6 


-76 


36 


15 


-2 


8 


19 


85 


3.5E-t-20 


2125 


Ab 





-56 


-38 






-44 


48 


-22 


48 




1.5E+21 



Note. — (1) KISSR ID. (2) Profile shape of Lya. Ab=absorption. P-Cygni=P-Cygni like. 2 Em=Two emission peaks. (3) Velocity offset between the 
centroid of the Lya line and Ha. (4)-(9) Velocity offset between the centroid of the intrinsic metal absorption line and Ha. (10) Mean velocity offset of 
intrinsic metal lines excluding Si III A1206. (11) Standard deviation of the velocities of the intrinsic metal lines excluding Si III A1206. (12) Velocity offset 
between Si ill A1206 and the mean of the lower ionization intrinsic metal lines. (13) Hi column density derived from Voigt profiles. 



- 25 - 



useful for comparing with global values derived from Lya line images, as they give an idea of the 
concentration of the emission and of how coincident it is with regions of high FUV surface bright- 
ness. Such Lya maps have been or will be obtained for a portion of our targets. We computed 
two values of the escape fraction, o ne using Ha and the other us ing H/3. Both adopt the case B 
recombination theory predictions of lOsterbrock &: Ferlandl ( 20061 ) for Tg ~ 10^ K and rig = 10^ 



cm~^, i.e., Lya/Ha=8.1 and Lya/H/3=23.1. L(Lya). The escape fractions were computed before 
and after correcting for the dust attenuation. We also computed the COS luminosity of the contin- 
uum at 1500 A, Li5oo, for comparison with the value from GALEX, which is given in Table [TJ In 
order to obtain -Lisoo; we extrapolated beyond the red edge of our data which is less than 1500 A. 
In this paper, L1500 is used for deriving the mass of the dominant stellar population and the star 
formation rate (see below). 

For correcting the hydrogen line fluxes and L1500 for the dust attenuation, we derived the total 
color excess, E{B — V)t, from the observed Ha/II/3 line intensity ratio, after correcting the Balmer 
lines for the underlying stellar absorption. The stellar abs orption correction w as done by adding 2 
A to the equivalent widths of the Balmer lines, following ( McCall et al.lll985l ). This increases the 
Ha flux by 2% and the H/3 flux by 13% (median values for the sample). We did not correct the 
interstellar Lya for the stellar contribution at Lya, since as mentioned in the introduction, this 
correction is expected to be small for stellar populations dominated by O and early B stars, as is 
the case for the Lya emitters in our sample. That O and early B stars dominate for the emitters 
is shown in section 3.5. The Lya emission line fluxes were obtained by numerical integration. The 
Balmer line fluxes were obtained from the SDSS data when available and from the KISS data 
otherwise. The SDSS spectra have higher spectral resolution. 

The total color excess is given by E{B - V)t = E{B - V)g + E{B - V)i, where E{B - V)g and 
E{B — V)i are the Galactic and i n trinsi c contributions, respectively. According to the Milky Way 
reddening maps of ISchlegel et al.l (Il998l ). E{B — V) q is small tow a rds ou r tar gets (see Table [H) 



thus most of the extinction must be intrinsic. Since 



Valls-GabaudI (|l993l ) and ICalzetti &: Kinnev 



(jl992l ) showed that by using a metallicity-dependent reddening law one could recover the theoretical 
Lya/H/3 ratio, we tried using a MW law for the spiral galaxies and an SMC law for the irregula r 
galaxies. For this purpose, we used the reddening laws tabulated in lOsterbrock &: Ferlandl ([200a). 
However, we foun d that the E(B — V\ values derived from these two laws are very similar. In 

fl), 



the end, and as in ICalzetti et al.l (|2000l ). we adopted a MW law for deriving E{B — V)t- For the 
theoretical Ha/H/3 ratio we used 2.86, which is the case B rec ombination theory prediction for 



Te = 10^ K and Ue = 100 cm s'^ of lOsterbrock Ferlandl (j2006|) 



In order to see if the relation E{B - V)s = 0.44 x E[B - V)t of ICalzetti et al.1 ([20001) holds 
for our data, we computed the color excess of the stellar continuum, E[B — V)s, usi ng the slope of 
the F UV continuum, /3. For obtaining E{B — V)s we used the method outlined in Iwofford et al.l 
(j201ll ). where the slope of the Milky- Way de-re ddened continuum is compared to that of a dust-free 
model of the appropriate metallicity. Following ICalzetti et al.l (120001 ). E{B — V)s was derived using 
the starburst reddening law. The comparison was done over the wavelength range from 1200 A to 
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the red edge of the data. The red edge is redshift-dependent, as show n in Figures 4 - 8, and it is 
always less than 1500 A. This range is different from what is used in ICalzetti et alj (2000), i-e., 
1250 - 1950 A. 

Our results are shown in Table [U which gives the KISSR ID of the galaxy (Column 1); the 
rest-frame value of EW(Lya) (Column 2); L(Lya), Lisoo from the COS spectra, L(Ha), and L(H/3) 
(Columns 3-6, uncorrected for dust extinction); the escape fractions within the COS aperture, i.e. 
Lya/Ha/8.1xl00 and Lya/H/3/23.1xl00 (Columns 7 & 8); E{B - V)t from the Balmer lines and 
E{B — V)s from the slope of the FUV continuum (Columns 9 &: 10); and the reddening-corrected 
values of L1500, Ha/H/3, and the escape fractions (Columns 11-14). 

In summary, we find that E{B — V)s = 1.27 x E{B — V)t (adopting the median of the sample), 
i.e., that the stellar continuum is m ore reddened than the nebular continuum, contrary to what 
was found by ICalzetti et al.l (120001). The differen t result is attributed to the larger optical and 



ultraviolet apertures used bv ICalzetti et al.l (|200Cl ). and to the different wavelength range used by 
the latter authors for fitting the UV continuum. We also find that the rest-frame EWs of the 
emitters are in the range 1 — 12 A, and that the reddening corrected escape fractions of the emitters 
are in the range 1 — 12%. The latter two results are further discussed in section 4. Finally, we found 
that L1500 from the COS data is lower by a factor of a few compared to the value from GALEX 
(compare column 4 of Table [6] with column 9 of Tabled]). This is attributed to the larger aperture 
of GALEX. 



3.4. SED Models 



We studied how the strength of Lya relates to the starburst phase and star formation rate 
of the dominant stellar population by comparing the COS observati ons with synthetic dust-free 



spectra computed with th e widely used package Starburst99 (S99, iLeitherer et al 



199S 



201c 



Vazquez Leithereijl2005l ) . The observed star formation rate is that within the COS aperture. The 
value of EW^Lya) depends on the meta llicity, IMF, star formation history, and evolutionary phase 



of th e stellar population (see models by IValls-Gabaudlll993l : ICharlot &: Falllll993l : IVerhamme et al 



20081 ). Beca use empirical spectr al stellar libraries in the UV are only available at the Galactic 
(Z^ = 0.013. ksphmd et al.ll2009l ) and LMC/SMC (Zlmc = 0.007 and Z^mc = 0.002. iMaeder et al 



I999I ) metallicities, we have to rely o n the oretical libraries at other metallicities. We used the the- 
oretical libraries of ILeitherer et al.l (I2OI0I ) corresponding to Z = 0.001, 0.004, 0.008, 0.020, and 
0.040. For comparison, we also computed models with the empirical Galactic and LMC/SMC stel- 
lar libraries. We will refer to simulations based on the theoretical and empirical libraries as the 
theoretical and empirical models, respectively. The theoretical models use stellar evolution tracks 
corresponding to Z = 0.020 and model atmospheres corresponding to Z = Zq. The empirical mod- 
els use tracks corresponding to Z = 0.020 or Z = 0.004. Our grid of models includes single stellar 
population (SSP) and continuous star formation (CSF) scenarios. Continuous star formation is 
more appropriate in our case, as the observed spectra include light from regions of a few kpc in size 



Table 6. Lya Escape Fractions 



Galaxy EW(Lya) logL(Lya) logLigoo logL(Ha) logL(H/J) /(Ha) /(H/3) E(B-V)t E(B-V), logLigoo Ha/H/3 /(Ha) /(H/3) 

A erg erg erg erg % % rnag mag erg % % 

uncorr uncorr uncorr uncorr uncorr uncorr Balmer /3 corr corr corr 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) (14) 



40 


< 




38.6 


40.02 


39.54 






0.04 


0.21 


39.5 


3.1 






108 


< 




39.1 


40.38 


39.85 






0.15 


0.17 


39.9 


3.5 






178 


< 




39.6 


41.21 


40.64 






0.23 


0.20 


40.5 


4 






182 


< 




39.1 


40.41 


39.88 






0.16 


0.19 


39.9 


3.6 






218 


< 




38.9 


40.49 


39.87 






0.33 


0.41 


40.6 


4.6 






242* 


10 


41.20 


40.1 


41.61 


41.07 


5 


6 


0.18 


0.21 


41.0 


3.7 


12 


14 


271* 


1 


38.74 


38.5 










0.17 


0.28 


39.7 








298* 


12 


40.16 


39.3 


40.65 


40.06 


6 


2 


0.28 


0.39 


41.0 


4.3 


12 


15 


326* 


4 


39.30 


38.7 


40.22 


39.61 


2 





0.16 


0.38 


40.3 


4.4 


7 


9 


1084* 


1 


38.99 


39.0 


41.03 


40.37 





1 


0.42 


0.44 


41.0 


5.2 


1 


1 


1567* 


3 


38.71 


38.3 


40.25 


39.63 





1 


0.15 


0.31 


39.6 


4.7 


2 


3 


1578* 


7 


41.00 


40.0 


41.55 


41.04 


3 


4 


0.12 


0.10 


40.4 


3.4 


6 


7 


1637 


< 




39.2 


40.43 


39.93 






0.10 


0.30 


40.5 


3.3 






1785 


< 




38.5 


40.38 


39.83 






0.18 


0.19 


39.3 


3.7 






1942 


< 




38.9 


40.00 


39.44 






0.21 


0.28 


40.1 


3.9 






2019 


< 




39.1 


40.63 


40.17 






0.01 


0.14 


39.7 


2.9 






2021 


< 




39.2 


40.59 


40.03 






0.21 


0.30 


40.5 


3.8 






2023 






38.6 


40.58 


39.91 






0.45 


0.31 


39.9 


5.5 






2110 


< 




39.1 


40.33 


.39.86 






0.00 


0.09 


39.5 


3 






2125 


< 




38.9 


40.73 


40.03 






0.49 


0.30 


40.2 


5.7 






Min 


1 


38.71 


38.3 


40.00 


39.44 








0.00 


0.09 


39.3 


2.9 


1 


1 


Max 


12 


41.20 


40.1 


41.61 


41.07 


6 


6 


0.49 


0.44 


41.0 


5.72 


12 


15 
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and such large regions must contain UV-bright star cl usters that sp an a range of ages. The SSP 
models are for comparison. We adopted a Kroupa IMF ( Kroupall200ll '). as at the high- mass end, it is 



considered to be universally applicable ( Bastian et al. 20ldl. We used the Geneva stellar evolution 



tracks for non-rotating single-stars with high-mass loss ([Scha ller et al.l Il992 ■ 



because models that account for stellar rotation and binarity ((Eld rid ge et al. 



Mevnet et al 



19941) 



20081 : iLevesque et al 



2OI2I ) are still under development and require calibration against observations over our range of 
metallicities. Finally, we did not compute models for times earlier than 2 Myr because this is the 
age of the youngest LMC main-sequence O stars, or greater than 30 Myr after the beginning of star 
formation since the SSP models are affected by incompleteness of the stellar libraries beyond this 
time and the continuum slope and age-sensitive stellar-wind lines of CSF models do not significantly 
change beyond this time. 



3.5. Starburst Phase, Stellar Mass, and Star Formation Rate 

The age of an FUV-bright stellar population can be derived spectroscopically by comparing 
the observed profiles of age sensitive lines originating in the winds of massive stars with model 
spectra. The strong Nv A1238.8, 1242.8. Si iv A1393.8. 1402.8 . and CiV A1548.2, 1550.8 resonant 



doub lets are widely used for this purpose (jTremonti et al.ll200ll : IChandar et al.ll2003l : IWofford et al 



20 111 ). In summary, the P-Cygni profiles of the Nv and the Civ doublets decrease in strength as 
the O stars evolve and expire. On the other hand, the Siiv doublet develops a P-Cygni profile 
in giant and supergiant O stars. Unfortunately, the Civ doublet and sometimes the Siiv doublet 
fall outside of the wavelength range that we observed with COS. In addition, the Nv doublet is 
sometimes contaminated with geocoronal emission. However, for targets showing clean and strong 
N y P-Cygni profiles, the age of the dominant population is younger than ~10 Myr (see figure 5 



m 



Wofford et al 



201lh. Fur t hermore, a large Ha equivalent also indicates the presence of a young 



population (lKennicuttill998l : iLeitherer et al.lll999l ). 



We compared the COS spectra with SSP and CSF models of 5 and 10 Myr. We select this 
age range because between 5 and 10 Myr, the strength of N v noticeably decreases and these ages 
are sufficient for our purpose of demonstrating how young the dominant stellar populations are. 
We found the best-fit model to the Nv and Siiv profiles by eye. The comparison was performed 
on the rectified spectra, as the models are dust-free but the data are not. The metallicity of the 
model is that most appropriate for the galaxy based on its oxygen abundanc^. Figures [T2] and [T3l 
show the comparisons of the observed N v and Si iv profiles with the theoretical models. Similarly, 
Figures [HI and [TSl show the comparisons with the empirical models. 

A comparison of the theoretical and empirical models in Figures [T2] and [T^ shows that at both 



^We adopted stellar evolution tracks with a metallicity loRmfZ/Zn ) ~ loKinf(0/H)/(0/H) 
corresponds to 12 + logj^Q(O/H)0 = 8.83, i.e., to the old solar value of lGrevesse fc Sauvall l|l99i 



where Zq = 0.020 
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ages and at all metallicities, Nv is stronger for the CSF models than for the SSP models, although 
the difference is less at the lowest metallicity (e.g., for 2019 and 2110). Furthermore, the theoretical 
CSF models yield stronger Nv profiles than the empirical models at similar ages and metallicities 
(e.g., 40 and 242). Conversely, the theoretical Siiv profiles show weaker emission components than 
the empirical ones. These differences do not severely affect our results. 

The observed N V profiles of Lya emitters 298, 326, and 1567 are contaminated with geocoro- 
nal Oi lines, and the Si IV profiles of 11 targets, including four of the seven emitters, are either 
incomplete or have low signal to noise. For this reason, and because the Siiv profiles from the 
theoretical models are too weak compared to observations at young ages (e.g., compared to the 
empirical models), in the rest of our age-dating analysis, we concentrate on the N V profiles. For 
targets with clean N v profiles, we find that the theoretical SSP models tend to be a better fit than 
the theoretical CSF models, except at the lowest metallicity (e.g., KISSR 2019 and 2110), where 
the SSP and CSF Nv profiles are very similar. This is true for the Lya emitters as well as for 
the Lya absorbers. One exception is 1084, which is better fitted by a CSF model. In any case, 
Lya absorbers such as spiral galaxy KISSR 218 and metal-poor galaxies KISSR 2019 and 2110 
have strong Nv profiles and large Ha EWs indicative of ages younger than ~10 Myr. Thus, they 
would be expected to have strong Lya emission. However, they are absorbers, which is consistent 
with their low gas-fiow velocities (see Table [5|). In addition, one cannot rule out attenuation of 
Lya due to the presence of dust. We also find that the emitters have Lya equivalent widths that 
are consistent with their young ages, but that are small compared to model predictions for young 
starbursts. Thus, the emitters also require the presence of scattering in Hi and attenuation by dust. 
The take away point is that even at young ages (< 10 Myr) targets can show Lya in absorption 
if the gas surrounding the starburst is almost static or if the dust column density surrounding the 
hot stars is significant. 

Table [7] gives the stellar masses and star formation rates derived from the UV spectra by 
adopting the ages in the second column of the table. For a given age, metallicity, IMF, and star 
formation history, the theoretical value of ^1500 only depends of the mass of the stellar population. 
We give two estimates of the mass and the star formation rate, one corresponding to the SSP 
and the other to the CSF model. The masses are derived from the comparison of the observed 
and computed luminosities of the stellar continuum at 1500 A (both dust-free). For the computed 
luminosity, we use the mean of the theoretical and empirical predictions at the adopted age and 
metallicity. Note that for the SSP models, it is common practice to define a characteristic star 
formation rate as the ratio of the stellar mass to the age of the population. This is what we quote 
in column 5 of Tabled For comparison, the table also gives the star formation rate d erived from the 



SPSS Ha luminosity (correc ted for reddening using the starburst attenuation law of lCalzetti et al 



(|2000l )) and equation (2) in iKennicuttI (|l998l ). We recall that the SDSS aperture is very similar 
to the COS aperture. This is not the case for the KISS data, therefore, we omit the optical SFR 
for targets without SDSS data. We also recall that the star formation rates are not global values. 
We find that the SSP models yield higher masses and star formation rates than the CSF models 
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by a factor of 1.4 (median). This is because CSF models predict larger values of L1500 than the 
SSP models at the adopted ages. We also find that the SSP models yield systematically higher star 
formation rates than those derived using Ha by a factor of 5 (median) with a standard deviation 
of 7. The UV SSP SFRs derived from L1500 should be viewed with caution as they are inversely 
proportional to the ages, which are expected to be underestimated by a factor of a few in cases 
where Nv is contaminated with geocoronal emission or does not show a strong P-Cygni profile. 
This is the case for KISSR 298 and 1637, which show the largest disagreement between the optical 
and the UV SFRs. In addition, the UV SSP SFRs depend on the masses, and Wofford et al. (2011) 
found systematically higher masses derived from Lisoo compared to masses derived from optical 
photometry for individual star clusters. On the other hand, in some cases, the agreement between 
the optical and UV SFRs is very good. This is the case for KISSR 1578, a target with a prominent 
Nv doublet and a high S/N spectrum. Unfortunately, only four emitters have Ha SFRs, which 
makes it difficult to study trends between Lya and the SFR. 



3.6. Wolf-Rayet Stars 



The detection of broad (FWHM ~ 10 A, Contll991 ) He li A4686 emission (hereafter, He il emis- 
sion) provides another way of estimating the age of a stellar population. In the absence of aii AGN , 
such emission originates in the dense winds of WR stars dBealsl I1929I: iBibbv fc Crowtheil 1201^ ) . 
whos e progenitors are rnassive stars (Minjtiai ^ 25, iMaeder &: Mevnetlll994l ') with lifetimes of ~ 5 
Myr (jMevnet k, Maedeilboosi ). Other optical WR emission hues are Cm A4650 and Niii A4640, 
which are generally weaker. For each galaxy, we show in Figure [16] the spectral region around the 
Hell A4686 line. The fi gure shows that in some c a ses, t he [Feiii] A4658 peak is taller than the 
Hell A4686 peak. As in iLopez-Sanchez &: EstebanI (j2010l ). we used multiple Gaussians to fit the 
spectral region in the vicinity of 4686 A. The Gaussian centers were fixed at wavelengths where 
lines were expected. We detect He 11 at the 2>a level or better in six cases, i.e., the amplitude of 
Hell is at least three times the standard deviation of the nearby continuum in six cases. In five of 
the latter cases, the He 11 line is broad, i.e., FWHM ~ 10 A. Unfortunately, for four of the seven 
Lya emitters, the optical spectrum is either unavailable (KISSR 271), of inadequate spectral reso- 
lution (KISSR 326 and KISSR 1567), or very noisy (KISSR 298). The three emitters with adequate 
optical spectra have He 11 detections with FWHM > 6 A. This includes KISSR 1084, which may be 
of composite spectral type, and whose broad He 11 emission is accompanied by narrow/nebular He 11 
emission. Note that for the latter galaxy, the He 11 peak is taller than the [Fe iii] A4658 peak, unlike 
in most galaxies with broad He 11 emission. On the other hand, two of the non-Lya emitters have 
He II detections and FWHM > 7 A. The lack of Lya emission in galaxies with broad He 11 emission 
could be due to the lack of significant gas flow in the direction towards the observer. This is the 
case of KISSR 178 and the highly i nclined spiral KISSR 2 1 25. F inally, the non-emitter KISSR 108 
only shows narrow Hell emission. IShirazi &: BrinchmannI (j2012l ) interpret the presence of narrow 
Hell emission without a broad component, as WR stars that are offset from the emitting gas. 
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1220 1240 1260 1220 1240 1260 1220 1240 1260 1220 1240 1260 1220 1240 1260 



Rest Wavelength (A) 

Fig. 12. — Comparison of observed and computed N v 1240 resonance doublets. The vertical gray 
lines give the rest wavelengths of the N v lines. The observations (thin black curves) are sometimes 
contaminated with geocoronal lines (dashed lines). We show SSP models (thick red curves) and 
CSF models (dotted-dashed blue curves) based on theoretical (T) stellar libraries. Such models 
only include stellar lines and are not optimized for the stellar Lya. The models correspond to 
either 5 or 10 Myr, depending on which model fits better. In each panel, we give the KISSR ID, 
the model age, and the model metallicity. The Lya emitters are marked with asterisks. 
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Fig. 13. — Similar to Figure [12] but for the Si iv 1400 resonance doublets. The vertical gray lines 
mark the positions of the Siiv lines. 
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Fig. 14. — Similar to Figure [12] but the models are based on empirical (E) stellar libraries. The 
empirical models are contaminated with Milky Way lines, in particular, Lya. Only two metallicities 
are available for the empirical models. 
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Fig. 15. — Similar to Figure [HI but for the Siiv 1400 resonance doublets. 
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Table 7. Stellar Population Properties 





Age 


A /T ^ nn 

iviass 


iviass 


or tx 


or tx 


or ix. 






OO-T 




oor 


Oor 






Myr 


IVl© 


iVi© 




iVi© yr 


J-vl© yr 


('-) 












K') 


40 


10 


7E+06 


3E+06 


0.7 


0.3 


0.1 


108 


5 


7E+06 


5E+06 


1.4 


0.9 


0.3 


178 


10 


7E+07 


3E+07 


6.8 


2.7 


2.6 


182 


10 


2E+07 


6E+06 


1.6 


0.6 


0.3 


218 


5 


4E+07 


2E+07 


7.2 


5.0 


0.7 


242* 


5 


9E+07 


6E+07 


17.6 


12.1 


5.6 


271* 


10 


lE+07 


5E+06 


1.1 


0.5 




298* 


10 


2E+08 


8E+07 


19.7 


7.9 


0.8 


326* 


5 


2E+07 


lE+07 


3.8 


2.6 




1084* 


5 


9E+07 


6E+07 


18.5 


12.8 


3.1 


1567* 


5 


4E+06 


3E+06 


0.8 


0.5 




1578* 


5 


3E+07 


2E+07 


5.0 


3.5 


4.0 


1637 


10 


6E+07 


3E+07 


6.5 


2.6 


0.3 


1785 


10 


5E+06 


2E+06 


0.5 


0.2 


0.3 


1942 


5 


lE+07 


8E+06 


2.3 


1.6 


0.2 


2019 


5 


5E+06 


3E+06 


1.0 


0.7 


0.3 


2021 


5 


3E+07 


2E+07 


5.2 


3.6 


0.6 


2023 


5 


7E+06 


5E+06 


1.5 


1.0 


1.2 


2110 


5 


3E+06 


2E+06 


0.6 


0.4 




2125 


5 


lE+07 


lE+07 


2.9 


2.0 


1.9 


Min 




3E+06 


2E+06 


0.5 


0.2 


0.1 


Max 




2E+08 


8E+07 


19.7 


12.8 


5.6 


Mean 




4E+07 


2E+07 


5.2 


3.1 


1.4 


Sigma 




5E+07 


2E+07 


6.2 


3.7 


1.6 



Note. — (1) KISSR ID identifier. The emitters are marked with asterisks. 
(2) Adopted age. (3) to (6) Stellar population masses and star formation 

rates corresponding to the COS aperture and derived by comparing the 
observed and computed UV spectra. The computed spectra correspond to 
the adopted ages aiid to SSP or CSF iiiodels. (7) Star formation rate derived 
from the reddening-corrected Ha luminosity. The last four rows give the 
minimum, maximum, mean, and standard deviation for each column. 
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4600 4650 4700 4750 4600 4650 4700 4750 4600 4650 4700 4750 4600 4650 4700 4750 4600 4650 4700 4750 



Rest Wavelength (A) 



Fig. 16. — Spectral region around Hell A4686. The legends give the KISSR ID and the ratios of 
the amplitude of Heii to the standard deviation of the nearby continuum. The spectra are from 
the SDSS DR7 except for galaxies KISSR 326, 1567, and 2110 whose spectra are from Salzer et al. 
(2005). We do not have an optical spectrum for 271. We degraded the SDSS spectral resolution 
to match that of the non-SDSS spectra, which is 2.4 A. The spectra are corrected for redshift and 
uncorrected for reddening. The dashed vertical lines mark the positions of the [Felll] A4658 and 
He II A4686 emission lines. We overlay in gray a fit to the continuum. 
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Table 8. Helium II 4686 A. 



Galaxy 


FWHM 


EW 


Comment 




4686, A 


4686, A 




(1) 


(2) 


(3) 


(4) 


40 








108 


2 


0.38±0.06 


narrow 


178 


7 


0.45±0.13 




182 








218 








242* 


6 


0.54±0.21 




271* 






no optical 


298* 








326* 






no SDSS 


1084* 


10 


1.19±0.04 


broad + narrow 


1567* 






no SDSS 


1578* 


9 


0.84±0.15 




1637 








1785 








1942 








2019 








2021 








2023 








2110 






no SDSS 


2125 


10 


1.26±0.23 




Min 


2 


0.38±0.06 




Max 


10 


1.26±0.23 





Note. — (1) KISSR ID. Wc mark the Lya emit- 
ters with asterisks. (2) & (3) FWHM and EW of Hen 
A4686 for galaxies with an SDSS spectrum and a 3o" 
detection. The last column indicates if the target has 
an optical spectrum or if a narrow He ii emission com- 
ponent is present. The last two rows give the minimun 
and maximum values in each column. 
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We compared the obs erved values of EW (Reii) with Starburst99 (jLeitherer et al.lll999l . l20ld : 



Vazquez &: Leithereijl2005l ) model predictions. Note that the observed values constitute lower limits 



due to the dilution of the Hen line in the stellar continuum of older stellar generations. Table [8] lists 
FWHM(He ii) and EW{}ie ii). Figure [T71 shows the comparison of observed and computed values of 
EW{}ieii) corresponding to single stellar populations (SSP) and continuous star formation (CSF) 
at the metallicities of the galaxies. Although Figure [17] shows that a given value of EW{Tieii) can 
be reached at multiple times, the range of times is restricted to 10^ — 10'' yr for the SSP models, 
and 10^ — 10^ yr for the CSF models. The presence of broad Heii emission rules out an IMF 
deficient in massive s tars or a popul a tion d ominated by late-B and A stars. Therefore, according 
to the predictions bv lCharlot &: Falll (jl993l ) one would expect positive values of EW{Lya) for the 
galaxies with broad He II emission. As previously mentioned, the lack of Lya emission in KISSR 
178 and KISSR 2125 is probably due to the dust geometry and inclination. In conclusion, we find 
no one to one correspondence between the Lya emission and the presence of WR stars. This is in 
agreement with our age-dating analysis based on the Nv profiles, where we concluded that some 
of the galaxies with young stellar populations are absorbers. 



4. Discussion 



4.1. Lya and Gas Flov^s 



Kunth et al.l (|l998l ) analyzed -fTST/GHRS observations of 8 nearby H ii galaxies covering metal- 
licities from 12-|-log(O/H)=8.0 to solar, and reddenings from E(B-V)=0.1 to 0.55. They found 
velocity offsets between the Hi and the Hii gas of up to 200 km s~^ in the four galaxies with Lya 
emission (P-Cygni profiles), whereas they found broad damped Lya absorptions at the velocity of 
the Hii gas in galaxies with no signs of Hi outflows from the LIS lines. These authors concluded 
that outflows are the main determining factor of the Lya escape from galaxies. Our kinematical 
results are in very good agreement with these findings. However, one still needs to assess the 
relative importance of scattering in H i and destruction by dust on the Lya escape. 

We also find that the Lya redshift is higher than the blueshift of the LIS line s by a factor of 
1.5 — 2. This is in agreement with the high redshift observations of IShaplev et al.l (|2003l ) and the 
expanding-shell radiation transport model of lVerhamme et al.l (|2006l ). 



Finally, note that the Lya troughs are the result of a superposition of Lya components. There- 
fore, the dominant componen t does not need to be at zero velocity. It is interesting that for our 
targets and for the targets in iKunth et al.l (|l998l ). the Lya absorption centroid is consistent with 
the velocity of the H ii gas. 
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Fig. 17. — Hell 4686 equivalent width versus time and metallicity for SSP (top) and CSF (bottom) 
models. In each panel, the dotted horizontal lines represent the upper and lower limits of EW(}ic ii) 
for the galaxies with broad Heii emission. We give the KISSR ID, SFH, and metallicity in each 
panel. The emitters are marked with asterisks. 
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4.2. KISSR 298 



KISSR 298 is the double-peaked Lya emitter. Such a profile is generally interpreted as a 
signat ure of Lyg radia,t i on tr ansfer through a non-static medium (see references given in section 
4.5 of IVerhamme et al.l ( 20081 ). For this object, we find an expansion velocity of 150 it 47 km s^^ 
from the LIS lines. This is high c ompared to the expansion velocities for the two double-peaked 



spectra in Verhamme et al. (j2008l ). which are almost consistent with a static medium (they are 



< 25 km s~^). 

In KISSR 298, the blueshifted peak is twice as strong as the redshift ed peak. Stronger blue 



peaks have been o bserved in star-forming galaxies at redshifts of z = 2 — 3 (IVerhamme et al 



2007 



Kulas et al.ll2012l ). They are qualitatively consistent with an infiow (jVerhamme et al.1 120061') and 
could be due to in-falling gas originating in a galactic fountain, as in the model of iTenorio-Tagle 
(|l99a ). The fact that in this system the dip between peaks essentially reaches the continuum is an 
interesting feature. Another possibility is that KISSR 298 has a broad Lya emission line blended 
with an absorption, and that the absorption center is shifted redward. The shift could be caused by 
the relative motion of the emission region and the absorb i ng ga s. A system corresponding to this 
description is shown in panel f of figure 3 of lWilman et al.l (|2005l ). which corresponds to the second 
cloud from the top in their figure 4. Unfortunately, the metal absorption lines of KISSR 298 are 
noisy and it is hard to tell if there is are unexpected absorptions near the wavelength corresponding 
to the Lya " absorption" . As can be seen in Figure [TTl the point equidistant to the two Lya peaks 
of KISSR 298, which is marked with a dashed line, has a velocity coincident with the Sin 1304 and 
C II A1336 absorptions. 



4.3. Lya versus Aperture Size and Distance 



We checked if there is an observational bias due to distance among the Lya emitters in our 
sample. We excluded KISSR 1084 and 1567 from the analysis as one is of composite type and for the 
other the TA failed. In order to separate the effect of distance and met allicity /reddening, we divided 
the emitters in two sets: the low metallicity irregulars (KISSR 242, 1578, and 271) and the higher 
metallicity spirals (KISSR 298 and 326). Each set spans a factor of about 2 in distance (1.75 for 
the spirals and 1.7 for the irregulars). The spirals have low inclinations. Therefore the inclination 
angle is a secondary effect. We find that L(1500), EW(Lya), and L{Lya) simultaneously increase 
with distance for both metallicity sets (see Table [6]). Note that this result is based on small number 
statistics and that the general trend at much higher redshifts is that brighter objects have lower 



equivalent widths (e.g. lAndo et al.ll2006l ) 



For the emitters, we find that EW(L ya) is low compared to values determined in other studies 
at similar redshifts with larger apertures. iGiavalisco et al.l ( 19961 ) found that EW{Lya) = 30 A on 
average for a sample of galaxies observed through the 10" x 2 0" aperture o f lUE . We attribute 
this to our smaller aperture. Indeed, as previously mentioned, lOstlin et al.l (j2009l ) found that in 
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nearby galaxies, the bulk of the Lya emission is in the diffuse component of the galaxies. We also 
find that the emi tters in our s ample have EW{Lya) slightly below the completeness limit of the 
z = 0.3 sample of lCowie et al.l ( 20ld ). which is 15 A. However, our sample has a median absolute B 
magnitude of -19.4 (see Tabled]), which is similar to that of the z = 0.3 sample (see their figure 23). 
This again could be an aperture effect. For Lyman-Break G alaxies (LBGs) at z ~ 3, the median 
value of EW{Lya) is zero (see figure 8 of lShapley et al.ll2003l ). Therefore, the EWs of the emitters 
in our sample are consistent with the latter median. 



4.4. Lyo versus Metallicity /Reddening 



If a homogeneous H I layer covers most of the H ii region, the value of EW(Lja) is expected to 
depend on the dust extinction because the Lya line flux is more strongly reduced (due to multiple 
scattering effects) than the adjacent continuum. 

Reddening due to the presence of dust is expected to increase with metallicity. We use oxygen 
as a gauge of metallicity. We re-deri ved the oxyg e n abu ndances of the 16 galaxies with SDSS 
spectra using the method described in lSalzer et al.l (120051 ). In summary, in cases where the weak 
[O III] A4363 line was not detected, we used the mean of the oxygen abundances derived from the 
EP and SDSS COARSE calibrations, and in cases where the R23 abundance could b e computed. 



we u sed the mean of the EP, SDSS and R23 abundances (all methods are defined in ISalzer et al 



20051 ). The oxygen abundances from the KISSR and the SDSS data are in agreement within the 
uncertainties and the final values are those given in Table [TJ 

We find that the two most metal-poor and least reddened galaxies, KISSR 2019 and 2110, 
show Lya in pure absorption, while the most metal-rich and most reddened galaxy, KISSR 298, 
has the highest value of EW{Lya). Strong Lya emission is expected from KISSR 2019 and 2110 
based on the strength of their Nv profiles and their high values of EW{Ha) (~ 500 A), which 
are characteristic of a young burst of star formation. However, these two targets show Lya in 
absorption. This is not surpr ising given that t his is also the case for I Zw 18, the most metal-poor 
galaxy in the local universe (lAtek et al.ll2009l ). Although we suffer from small number statistics, 
the lack of one to one corre spondence between the strength of Lya and metallicity /reddening is in 
agreement with the work of Giavalisco et al. (jl996l ) , which is based on larger optical and ultraviolet 
apertures than what we use in the present study. Therefore, it is not clear that larger aperture 
measurements yield a different result regarding Lya versus metallicity. Our upcoming Lya plus 
Ha /H/3 ratio maps of three of the galaxies in our sample should help establish the role of dust and 
metallicity in regulating the Lya escape. These maps will be obtained as part of HST GO-12951 
and will enable the measurement of spatially resolved and integrated Lya and reddening values. 
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5. SUMMARY AND CONCLUSION 

1. We used the HST COS to observe 20 nearby (z ~ 0.03) galaxies spanning a broad range 
in properties (Tabled]). Our data covers the wavelength range 1138 — 1457 A and samples 
circular regions of 0.9-2.4 kpc in radius within each galaxy. We studied correlations between 
the Lya line properties and galaxy properties including distance, metallicity, reddening, star- 
burst phase, and gas flow properties. Our study includes the analysis of ancillary optical 
spectroscopy from SDSS and KISS. Based on the BPT diagram (Figured]), the galaxies are 
dominated by star-formation except for KISSR 1084, which is of composite spectral type. 

2. We found Lya pro file shapes representative of what has been observed betwe en redshifts of z ~ 



and z ~ 5 (e.g. jKunth et al.lll998l : IXapken et al.l 120071 : iKulas et al.ll2012l ). including strong 



absorptions, P-Cygni profiles with redshifted emission, and double emission (Figure [TU]) . 

3. Seven out of the 20 galaxies have net Lya emission, including spiral and irregular galaxies and 
composite galaxy KISSR 1084 (Figures [3] andfTO]). For these seven emitters, the Lya equivalent 
width is in the range 1 - 12 A and the range of escape fractions within the COS aperture is 
1 — 12 % (Table [6]). The comparison with escape fractions at higher redshifts is complicated 
by the fact that at higher redshifts, observations include more if not all of the galaxy in the 
aperture. 

4. A face-on spiral, KISSR 298, shows two peaks of Lya emission separated by 370 km s~^, one 
blueshifted and one redshifted and half as strong. Unfortunately, we were unable to determine 
if this profile is a strong emission with absorption within, as LIS lines for this galaxy are too 
noisy. 

5. Excluding KISSR 298, the emitters have Lya peaks redshifted to 172 it 49 km s~^ with respect 
to Ha (mean it standard deviation). Including all emitters, the mean O i A1302 and C ii A1334 
absorption blueshifts are -117 ± 47 and -94 ± 20 km s'^ Therefore, excluding KISSR 298, 
the Lya redshift is larger than the interstellar absorption blueshift by approximately a factor 
of 1.5 — 2, as found at high redshift (e.g. Shapley et al. 2003), and in agreement with radiation 
transfer models (e.g., Verhamme et al. 2006). 

6. For the absorbers, Lya is centered at 0ib38 km s^^ and O i and C ii have centroids at —23 it 28 
and —21 it 31 km s~^. Thus, most absorbers are consistent with having static or low velocity 
Hi gas. This supports earlier findings (e.g. of Kunth et al. 1998). 

7. The outflow velocity of Si iii A1206 is typically, 30 — 50 km s~^ greater than that of the H i gas 
(Tabled]). 

8. We found no one to one correspondence between the strength of Lya and the redden- 
ing/metallicity measured within the COS aperture (Tables [6] and [7]) . In particular, we found 
three low- inclination spirals each showing different Lya profiles, i.e., pure absorption, double 
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emission, and P-Cygni profile, and the two most metal-poor and least-reddened galaxies in 
our sample are pure absorbers. 

9. We find stellar populations consistent with SSPs of ages in the range 5-10 Myr among the 
Lja emitters and the Lya absorbers. Similarly, we detect signatures of the presence of Wolf- 
Rayet stars among the emitters and the absorbers. Therefore, a young starburst does not 
guarantee strong Lya emission. 

10. In conclusion, a picture emerges where the Lya photons escape through regions of low H I and 
dust column densities, where gas outflows are occurring. 

11. In HST cycle 20, we will map three low- inclination spirals from our sample in the Lya line 
and the H/3/Ha ratio (proxy for dust attenuation). This will make it possible to study the 
relative importance of gas flows and dust content in regulating the Lya escape. 
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